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Abstract 

This paper describes a joint activity involving 
NASA and Army researchers at the NASA Langley 
Research Center to develop optimization procedures 
to improve the rotor blade design process by integrat- 
ing appropriate disciplines and accounting for all of 
the important interactions among die disciplines. 
The disciplines involved include rotor aerodynamics, 
rotor dynamics, rotor structures, airframe dynamics, 
and acoustics. The work is focused on combining 
these five key disciplines in an optimization proce- 
dure capable of designing a rotor system to satisfy 
multidisciplinary design requirements. Fundamental 
to the plan is a three-phased approach. In phase 1, 
the disciplines of blade dynamics, blade aerodynamics, 
and blade structure are closely coupled while acoustics 
and airframe dynamics are decoupled and are accounted 
for as effective constraints on the design for the first 
three disciplines. In phase 2, acoustics is integrated 
with the first three disciplines. Finally, in phase 3, 
airframe dynamics is integrated with the other lour 
disciplines. Representative results from work per- 
formed to date are described, lhese include optimal 
placement of tuning masses for reduction of blade vi- 
bratory shear forces, integrated aerodynamic/dynamic 
optimization, and integrated acrodyniimi c/dynam- 
ic/structural optimization. Examples of validating the 
procedures are described. 


* Deputy Head, Interdisciplimiry Rese;irch Office, 
Associate Fellow AIAA, Member AltS, ASMF 

**Researeh Engineer, Interdisciplinary Research 
Office, Member AIAA, A1 IS 
t Research Engineer, AVSOOM Acrostruclures 
Directorate, Member AHS 


Introduction 

An emerging trend in the analytical design of air- 
craft is the integration of all appropriate disciplines in 
the design process. 1 ' 2 This means not only 
including limitations on the design from the various 
disciplines, but also defining and accounting for 
interactions so that the disciplines influence design 
decisions simultaneously rather than sequentially. In 
rotorcrafl design, the appropriate disciplines include 
aerodynamics, dynamics, structures and acoustics. 
This paper describes an activity for developing the 
logic elements for helicopter rotor design optimiza- 
tion which includes the above disciplines in an inte- 
grated manner. 

Rotorcraft design is an ideal application for inte- 
grated multidisciplinary optimization. In current de- 
sign practice however, the process has generally been 
sequential, (i.e. single-discipline oriented) 3 ' 6 rather 
Uuui integrated. In early 1985, several occurrences led 
the NASA Langley Research Center to address the 
multidisciplinary design problem. The Inter- 
disciplinary Research Office within the NASA 
I .anglcy Structures Directorate was charged with the 
development of integrated multidisciplinary optimiza- 
tion methods. Nearly concurrently, the Army 
Acrosiructurcs Directorate at Langley established the 
goal of improving rotorcraft design methodology by 
"discipline integration." Close cooperation between 
the NASA and Army organizations led to plans for a 
comprehensive research program for an integrated ana- 
lytical design capability. As a result of the common 
goals in rotorcraft design, a group of NASA/Army re- 
scarchcrs (referred to herein as die Langley Team) ini- 
tiated the approach described in this paper. 

The focus of the NASA/Army research is to 
develop strategics, logic, and formulations for 
integrated multidisciplinary design optimization will) 
a view toward their application in the rotorcraft 
community. Toward this end, the Langley Team 
produced a comprehensive plan 7 which was sent to 
each of the principal helicopter companies in die 
United Stales to obtain their critique. Subsequently, 
the I anglcy Team completed a lour of the companies 



to review the results of the critiques. An updated plan 
was produced which was more representative of 
industry design philosophies and guides the ongoing 
research. The purpose of this article is to describe the 
plan and to summarize the most recent results and 
progress in the work. Specifically, the paper contains 
a summary of the plan, the sequence of the develop- 
ment, and some recent results. These results include 
optimal placement of tuning masses for reduction of 
blade vibratory shear forces, integrated aerody- 
namic/dynamic optimization, and formulation for in- 
tegrated aerodynamic/dynamic/structural optimization. 

General Approach a nd Scope 

Development Strategy 

A three-phased approach for the activity is illus- 
trated in figure 1. In phase 1 the rotor blade 
aerodynamics, dynamics and structural analyses are 
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Fig. 1 Three-phase development strategy 
for rotorcraft optimization 

coupled and driven by the optimizer. The 
optimization of the blade aerodynamic geometry, 
blade stiffness and mass distributions and the detailed 
internal structure is performed. The design 
requirements influenced by the acoustics and airframe 
dynamics are accounted for indirectly by constraints 
on the blade design. In phase 2 the acoustic analysis 
is fully integrated with the blade aerodynamics, 
dynamics, and structural analysis. Constraints will 
include limits on sound pressure levels and 
directivity. The design produced in phase 2 will 
satisfy acoustics goals. Airframe dynamics in phase 
2, as in phase 1, is accounted for by effective 
constraints on the blade dynamics, aerodynamics and 
structural behavior. Finally, in phase 3, airframe 
dynamics is integrated and the result is a fully 
integrated optimization strategy. 

Sequence of Tasks for l’hase 1 

Figure 2 depicts the general sequence of single 
and dual discipline optimization tasks diat lead to a 


lully-integrated rotor blade aerodynamic/dynam- 
ic/structural optimization procedure. The dynamic 
optimization bubble in the figure represents the work 
from references 5, 6, and 8- 1 1 . The aerodynamic per- 



Fig. 2 Sequence of optimization elements 
for phase 1 


lbrmancc optimization bubble is based on the work 
described in reference 1 2. The left-center bubble 
represents the integration of an aerodynamic loads 
analysis with dynamics; a procedure wherein the 
airloads can be adjusted by changes in the design 
variables to reduce dynamic response. This work is 
described in reference 13. The arrows leading from 
these three bubbles indicate that a merger of these 
three produced a fully integrated aerodynamic/dynamic 
procedure described in reference 14. The structural 
optimization bubble indicates that some work in rotor 
structural optimization was being initiated 
simultaneously with the integrated work and is 
discussed in reference 15. The strategy for integration 
of the structural optimization with the dynamics and 
aerodynamics represented by the bubble on the far 
right is a multilevel formulation based on the theory 
described in references 16 and 17. In the present for- 
mulation, the planform shape, pretwist, stiffnesses 
and mass distributions of the blade are determined in 
the upper level and the detailed sizing of the blade 
spar takes place in the lower level. 


Status/Results 


1'he development of the strategies for the ele- 
ments in figure 2 is nearly complete. Validated 
demonstrations of optimum placement of tuning 
masses for vibration reduction, integrated aerody- 
namic/dynamic optimization, and the strategy for in- 
tegrated aerodynamic/dynamic/structural optimization 
;tre highlighted in this section of the paper. 
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Optimum Location s of Vifrf-»«ion Tuning 

Masses 

The objective of this work was to develop and 
validate a dynamic optimization procedure that 
systematically determines the best values and loca- 
tions for tuning masses for reducing vibratory vertical 
hub shear. Optimal placement of the mass tailors the 
mode shapes and airloads thus reducing generalized 
force and response of the blade. The method entails 
formulating an optimization procedure that employs 
the tuning masses and their locations as design vari- 
ables to minimize vertical hub shear with the small- 
est possible mass penalty. Equation (1) defines the 
objective function, f, which is a combination of 
vertical hub shear and added mass where K is the 
number of shear harmonics to be included in the 
objective function, Mj is the jth tuning mass, and 
Sref * s a reference value of vertical hub shear. 
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The additional design variables, Pk * which ap- 
pear in the objective function (equation 1) and the 
constraints shown below in equation (2) are "pseudo 
upper limits" on the calculated shear harmonic 
amplitudes, S^, 


As described in detail in reference 10, the im- 
plementation of the approach was to combine the op- 
timizer CONMIN 18 with CAMRAD/JA 19 . The 
latter code was used to calculate mode shapes, 
frequencies, airloads and hub shears. The use of 
CAMRAD/JA enables the variation in the airloads 
due to changes in the design variables to be taken into 
account 

The bar chart in figure 3 compares the results of 
applying the optimization procedure to place three 
masses and six masses on a one-sixth, Mach-scaled 


COMPARISON OF BASELINE AND OPTIMIZED 
DESIGNS FOR 3-MASS AND 6-MASS CASES 
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Fig. 3 Results for optimal placement and 
sizing of three and six masses 


S k /p k -lS0 k = 1,2,. ...K 
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These constraints express the requirement that S k be 
less than the value of p k Consequently, the opti- 
mizer will tend to increase the values of p k to satisfy 
the constraints but at the same time will attempt to 
decrease the values of P k to minimize the objective 
function. This results in a compromise on the values 
of Pk which forces a reduction in the values of S k 
thus reducing the hub shear harmonics while incur- 
ring the smallest possible mass penalty- Additional 
constraints include upper and lower bounds on the 
first two flapwise frequencies of the blade to avoid 
resonance as shown in equation (3) 


(Oj /(O ui -l<0| 

1 - Wj / o>|j < 0 J 


i = 1,2, 


...,l 
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where oij = (Oj / Q, (Oj is the ith natural frequency, 
o) U i and (0|j are the upper and lower bounds on the 
frequency, respectively, (1 is the blade rotational 
speed, 1 is the number of constrained frequencies. 


model of the growth version of the U. S. Army's UH- 
60A helicopter rotor blade. The axis on the left 
measures the shear harmonic amplitudes and the axis 
on the right measures the amount of tuning mass that 
was added. The groups of columns from left to right 
represent the third, fourth, and fifth harmonics of 
shear and the mass, respectively. The columns 
within each group from left to right represent the 
values of the baseline which has no added mass, the 
three-mass case and the six-mass case, respectively. 
Using three masses, the procedure was able to reduce 
the shear by about eight percent for the third and 
fourth harmonics and four percent for the fifth with a 
ten percent increase in total blade weight. Using six 
masses, the third harmonic was reduced by 24 percent 
from the nominal, the fourth harmonic was reduced 
34 percent, and the fifth harmonic 32 percent. The 
amount of mass that was needed to achieve these 
reductions was approximately a 30 percent increase in 
the total blade weight which would probably 
disqualify this blade in practice. However, this result 
shows that the method trades shear reduction for mass 
and it confirms the hypothesis that large reductions in 
vibratory shear require a large mass penalty. 
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As described in reference 10, the procedure has 
also been applied to a blade test article (shown in fig- 
ure 4) that has the capability for adding tuning masses 



Fig. 4 Test article for experimental vali- 
dation of mass placement method 


predictions were respectively eleven percent and 
twelve percent below the range. This is fairly good 
agreement considering the difficulty of predicting 
fixed-systems loads with existing analysis codes. 

Integrated Aerodvn a mic/dvnamlf 
Optimizati on fAID 

A fully integrated aerodynamic/dynamic optimiza- 
tion procedure 14 which combines performance and 
dynamics analyses for hover, forward flight, and 
maneuver with a general purpose optimizer has been 
developed. The maneuver flight condition simulates a 
sustained pull-up maneuver in terms of a load factor 
on die forward flight lift requirement. 

O b j ect i ve — Function - The procedure mini- 
mizes an objective function which is a composite 
measure of performance and dynamics. Specifically, 
die objective function is a linear combination of 
power required (for hover, forward flight, and maneu- 
ver) and vibratory hub shear 


along the blade span. The test article, designed to 
study passive means for minimizing fixed-system 
loads, was tested in the Langley Transonic Dynamics 
Tunnel. The test data included 4/rev hub shear as a 
function of the location of a single tuning mass of 
0.27 lbm for several flight conditions. The graph in 
figure 5 shows the comparison between the 
optimization results and the test data for three advance 
ratios (0.25, 0.30, and 0.35). The test data is shown 
as a band since the data was available only at ten per- 
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where IIP^, HP^, and HP m are the powers required in 
hover, forward flight, and maneuver, respectively. 
S N ff is a,c N P 4 * rev nonrotating vertical hub shear 

in forward flight and N is the number of blades. The 
terms k j , k 2 > k^, and k^ are weighting factors. 

"'’•‘ref’ HPff ref’ HPm re f’ aDd S N ref reference 
values used to normalize and nondimensionalize the 
objective function components. 

Design Variables . The design variables are 
shown in figure 6 and consist of aerodynamic quanti- 
ties describing the blade planform and pretwist and of 
dynamic quantities describing the blade structural 
properties. The four aerodynamic design variables are 


Fig. 5 Comparison of optimum locations 
of single mass with test data 

cent increments along the span. For the 0.35 advance 
ratio case, the optimization procedure predicted an 
optimum location that was within the range of die 
test data. For the other two cases ( 0.25 and 0.30 ) the 
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Fig. 6 Design model for integrated aero- 
dynamic/dynamic optimization 

the point of taper initiation y (r root chord c r , taper 
ratio Cj/Cj, and maximum prelwist 0 lw . Tlie blade is 
rectangular to and then tapers linearly to the tip. 

The pretwist (blade structural and aerodynamic twist 
are assumed to be the same) varies linearly from the 
center of rotation to the tip. Nine dynamic design 
variables include the blade stiffnesses in Hie lagwise 
(inplane) and flapwise (vertical) directions (denoted by 
El/ and Elf) and torsional stiffnesses (denoted by GJ) 
at the blade root, point of taper initiation, and blade 
tip. The stiffnesses are assumed to vary linearly 
between these three points. Ihe remaining six 
dynamic design variables are three tuning masses 
(denoted by m j, m2, and m3) and their distances from 
the center of rotation (denoted by yj, y2. and y3). 
The total blade mass consists of the structural mass 
(which remains constant) plus the sum of the tuning 
masses. 

Constraints • The constraints arc grouped into 
performance constraints and dynamic constraints. The 
performance constraints are imposed for all three 
flight conditions. The dynamic constraints are im- 
posed only in forward flight and maneuver. Hie blade 
is designed for a constant lift in forward flight and a 
constant lift in maneuver. 

The performance constraints are on power re- 
quired. stall, trim, and blade tip chord. The limit on 
power is that the power required in hover, forward 
flight, and maneuver be less than the power available. 
The requirements that the airfoil sections not stall and 
that the drag divergence Mach number arc avoided are 
expressed as upper limit constraints on lire airfoil sec- 
tion drag coefficients as functions of the angle of at- 
tack and Mach number. These constraints arc evalu- 
ated at every 15 degrees around the azimuth in forward 


flight and maneuver. An isolated rotor analysis is 
used which trims the rotor to constant lift and drag 
and zero flapping angle relative to the shaft using col- 
lective, lateral cyclic and longitudinal cyclic pitch. 
Trimming to a constant lift ensures that the rotor re- 
tains tire specified lift capability even if solidity de- 
creases. The trim constraint is implemented in for- 
ward flight and maneuver. The final performance 
requirement is a lower limit on the blade tip chord. 

The dynamic constraints are “windows” on fre- 
quencies, total blade weight, and autorotalional iner- 
tia. The constraint on frequency (either a bending or a 
torsional frequency) avoids integer multiples of the 
rotor speed £ 2 . The constraint on blade mass sets the 
maximum value (the total blade mass is the sum of 
die constant structural mass and the tuning masses) 
and die constraint on autorotational inertia requires 
safe autorotation in case of engine failure. 

Rotor Analyses - The analyses used in this 
work are Hie Langley-developed hover analysis pro- 
gram HOVT (a strip theory momentum analysis 
based on reference 20) and the comprehensive 
helicopter analysis program CAMRAD/JA for 
forward flight and maneuver. HOVT is used to 
predict power required in hover using nonuniform 
inflow (no wake is included). CAMRAD/JA is used 
to calculate rotor performance, loads, and frequencies. 
In this work the CAMRAD/JA analyses are 
performed with uniform inflow with empirical inflow 
correction factors. Both HOVT and CAMRAD/JA 
use tables of experimental two-dimensional airfoil 
data. 

Optimization Methods - The codes used for 
optimization arc the general purpose optimization 
program CONMIN and an approximate analysis used 
to reduce the number of HOVT and CAMRAD/JA 
analyses during the iteration process. CONMIN is a 
general purpose optimization program which uses the 
method of usable-feasible directions for constrained 
function minimization. The approximate analysis is 
used to extrapolate the objective function and 
constraints with linear Taylor Series expansions 
using derivatives of the objective function and 
constraints with respect to die design variables. The 
assumption of linearity is valid over a suitably small 
change in the design variable values and will not 
introduce a large error into the analysis provided the 
cltangcs are small. Errors which may be introduced 
by use of the approximate analysis are controlled by 
imposing “move limits” on each design variable 
during die iteration process. A move limit which is 
specified as a fractional change of each design variable 
value is imposed as an upper and lower design 
variable bound. 
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Im plementation of Optimisatio n Procedure . 
The optimization procedure (figure 7) consists of an 
outer loop denoted by “Cycle” and an inner loop de- 
noted by “Iteration”. First, preassigned parameters 
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Fig. 7 Flowchart for integrated aerody- 
namic-dynamic optimization 

such as the blade radius, airfoil distribution, and 
number of blades are set. An optimization cycle is 
initiated. The aerodynamic and structural properties 
such as twist and chord distributions, radial station 
locations, solidity, blade weight, and autorotational 
inertia are calculated using the current design variable 
values in the box labelled “Design variable 
preprocessors”. The HOVT analysis is then per- 
formed to obtain the power required in hover. Two 
CAMRAD/JA analyses (forward flight and maneuver) 
are then performed to obtain the power required, trim 
information, coefficients of drag for the stall con- 
straints, natural frequencies, and hub shears. This in- 
formation is then used to formulate the objective 
function and constraints. Since CONMIN and die 
approximate analysis need derivatives of the objective 
function and constraints, a sensitivity analysis is per- 
formed to obtain finite difference derivatives of the 
objective function and constraints with respect to the 
design variables. These derivatives are obtained by 
perturbing each design variable one at a time and go- 
ing through the design variable preprocessor, HOVT, 
and CAMRAD/IA analyses. The inner loop consists 
of CONMIN and the approximate analysis. New val- 
ues for the design variables are obtained and the outer 
loop is re-entered. Convergence is obtained if the ob- 
jective functions from three consecutive cycles are the 

same within a tolerance of 0.5 x I O' 5 . 

Va l idation by Comparison with 

Conventionally-designed ki.h. 

- The proce- 
dure was tested from an arbitrary starling point to 
determine how well it could reproduce the design of 
an existing wind tunnel model of a growth utility 


rotor blade which was designed previously by 
conventional methods (not using formal optimization 
techniques). This existing blade will be referred to as 
the actual blade. This blade has a rectangular plan- 
form to 0.80R (80 percent radius) and then tapers to 
the up with a 3-to-l taper ratio. The blade has a ra- 
dius of 56.22 inches and a root chord of 5.40 inches 


I lie flight conditions are (1) a constant lift of 
scaled l-g (331 pounds), propulsive force of 32 
pounds, and an advance ratio of 0.35 for the forward 
flight condition and (2) a constant lift of 401 pounds 
a propulsive force of 23 pounds, and an advance ratio 
or 0.3 for the maneuver flight condition. The 
maneuver flight condition has a load factor of 1 22 
Since f, corresponds to a rigid body mode and f 2 is 
the I per rev, the first two frequencies are not 
constrained. Constraints are placed on the first four 
bending frequencies (f 3 through f 6 ) and the first two 
torsional frequencies (t, representing the rigid body 
torsional mode due to the control system stiffness and 
t 2 representing the first elastic torsional mode). The 
analysis of the blade calculates a total weight of 3.05 
lbs and an autorotational inertia value of 341 1 lbm- 
m The blade is to be designed so that the weight is 
not increased by more than 15 percent and the 
autorotational inertia is increased by at least 1 percent 
from that of the actual blade. The values for 
minimum tip chord, power available, and maximum 
allowable drag coefficient are 1 in, 20 hp, and 0.1, 
respectively. A frequency window of ±0. 1 per rev is 
used for the frequency constraints. 


The initial trial blade design (the starting point for 
the optimization) is a blade which has a rectangular 
pianfonn with a maximum pretwist of -9.0 degrees 
and blade root chord of 5.40 in. This blade has the 
same root chord, mass distribution, and stiffness dis- 
tributions at the root, 0.8R (point of taper initiation 
of the actual blade), and the tip as the actual blade. 

1 he stiffnesses are assumed to vary linearly between 
tiiesc points. The nonstructural mass distribution de- 
pends only on the tuning masses and their locations. 
Note dial the initial blade does not satisfy the mini- 
mum autorotational inertia requirement. 

Since a four-bladed rotor is used as the test prob- 
lem, die 4 per rev nonrotating hub shear is used for 
S N (| in ,llc objective function given by equation 4. 

Ihc normalizing factors HP h , HP ff , and HP 

n ref n rcf m ref 

are each chosen to be 15 hp and Sm is chosen to be 

'ref 

2 Ibf (based on analysis of the initial blade). The 
actual blade was originally designed for performance. 
Therefore, the objective function is chosen to be one 
dominated by performance with little emphasis on 
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dynamics. Of the three flight conditions, it is as- 
sumed that it is most important to reduce the power 
required in hover - therefore, this term will have twice 
the weight as the other two horsepower terms. 
Several values were tried for the weighting factor on 
the hub shear term. It was found that to obtain the 
proper balance between performance and dynamics, IC 4 
has to be about three orders of magnitude smaller than 
kj. Thus, for this case, the weighting factors arc 
chosen to be kj =15.0, k 2 =k 3 = 7 . 5 , and k 4 = 0.025. 


Figure 8 includes the shape and response for the 
initial trial, final, and actual blades. The Figure 
shows comparison of performance and dynamics mea- 

InlUal trial Final Actual blada 

design design 

% I ■ v. 7 


Hover hp 

15.35 

14.41 

14.64 

Forward flight hp 

13.52 

12.54 

13.13 

Manauvar hp 

12.26 

11.78 

11.83 

4 par rev 
hub ahaar Ibf 

2.19 

1.17 

1.52 


Fig. 8 Comparison of optimization-based 
design with actual blade 

sures (horsepower for hover, forward flight, and 
maneuver and the 4 per rev forward flight vertical hub 
shear). The final design has the same pretwist as the 
actual blade. Both planfomis are similar with die 
final design having less solidity than the actual blade. 
Specifically, the root chord was 4.4 in, the taper ratio 
was 1 . 8 , the point of taper initiation was 0 . 68 , and 
the maximum pretwist was -16 degrees. For die 
actual blade these values were 5.4 in, 3.0, 0.80, and 
-16 degrees, respectively. The difference in planforms 
is primarily due to the choice of flight conditions. 
The actual blade was designed by use of parametric 
studies for slightly different flight conditions and 
design requirements. 

Ya liditivn By U st p f Test A rtick 

The AD procedure has been used to design a lest 
article which will be used to validate die procedure 
through wind-tunnel testing. Along with a baseline 
blade, the validation blade will be tested for aerody- 
namic and dynamic performance in the Iran sonic 
Dynamics Tunnel (TOT) at Langley. Ibe baseline 
blade was previously designed by the U.S. Army 
Aerostructures Directorate at die Langley Research 
Center and was made available to die Langley learn 
for use in the current work. The purpose of the tests 
will be twofold: to verify dial die optimization proce- 
dure can produce a blade that will improve the aerody- 
namic and dynamic performance relative to a baseline 


design; and to verify the behavior and trends predicted 
by die analytical procedures. 

The validation blade was designed for a five- 
bladed rotor for the following conditions: in hover, a 
constant 1 .0 g scaled lift of 301 lb; in vertical rate of 
climb of 1000 ft/min; in forward flight, a constant 
1.0 g lift of 301 lb and a constant propulsive force of 
28 lbs at an advance ratio of 0.326; a maneuver simu- 
lated by a load factor of 2.0 g (602 lb) at an advance 
ratio of 0.233 with a propulsive force of 9.0 lbs. The 
objective function is that shown in equation ( 4 ) 
with the following weighting factors: kj = 1.0 kj 

= 1.0 k^ = 0.5 k^ = 1.0. The optimization proce- 
dure was carried out with the initial blade design be- 
ing the baseline blade which has a rectangular plan- 
fomi and pretwist of -8.0 degrees. The CAMRAD/JA 
analysis employed in die analyses used a wake model 
with rigid geometry with vertical convection by the 
mean inflow. 
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Fig. 9 Validation test article 

Figure 9 shows die results of the optimization. 
The analytical predictions in the upper Table show a 
predicted improvement in hover power, forward flight 
power and maneuver power of six percent, three per- 
cent, and four percent, respectively. The improvement 
in vibratory load level is predicted to be minimal. 
The validation last article is presendy in the fabrica- 
tion stage and tesdng is expected to be initiated dur- 
ing calendar year 1993. 

Integrated Aerodvnamic-Dynamic-Structural 
O ptimization (ADS) 

Multilevel Optimization Strategy 

The ADS optimization strategy is based on the 
method of multilevel decomposition described in ref- 
erences 16 and 17. In this case there are two levels. 
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In the upper level (figure 10) the goal is to optimize a 
combination of aerodynamic and dynamic performance 
while satisfying constraints on the aerodynamic, dy- 
namic and global structural behavior (i.e., average 
strain). The upper level is essentially the same 
as the AD 



Fig. 10 Upper level of integrated aerody- 
namic-dynamic-structural optimization 


design variables from the upper level are passed to the 
lower level. The lower level objective function for 
each section is a measure of how close the El's, GJ 
and M are to those required at the upper level. 
Constraints are enforced at each section to assure that 
the lower level produces an EA which is at least as 
large as that required in the upper level. Constraints 
arc also enforced on the axial and shear stresses in the 
blade.Thc design variables are the detailed dimensions 
of the blade cross sections. 

Example of ADS Procedure 

An example of the ADS optimization procedure 
is described in this section of the paper. For 
illustrative purposes and simplicity, it is assumed 
that the blade is constructed from an isotropic 
material. Although no numerical results from this 
formulation are included in the paper, this example is 
the basis for numerical calculations which are in 
progress for demonstrating the methodology. 


procedure with the addition of strain constraints, 
design variables representing section masses M 
(including structural and non-structural mass but 
excluding tuning masses), extensional stiffness 
variables EA at several spanwise locations, and 
coordination constraints which link the upper and 
lower levels. 

In the lower level (figure 11), the sizing of the 
internal blade structure takes place. The purpose of 
the lower level optimization is to assure that a 
structure can be sized to provide the required 
stiffnesses and section masses needed tit the upper 
level and also to assure the structural integrity of the 
blade. The lower level optimization is performed in 
parallel for several spanwise cross sections. As 
indicated in figure 11, the HI*. GJ*, EA*, and M* 
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Upper Level Obje ctive Function - The 
objective function is a combination of rotor horse- 
power and transmitted vibratory loads at several flight 
conditions 
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where F is the objective function, k| are weighting 
factors, HPj are horsepowers for three flight condi- 
tions listed below, is the N/rev vertical hub 

shear in forward flight, and HP ref andS N arerefer- 

1 ref 

ence values of horsepower and load respectively. The 
set of flight conditions is: 



Flight 

condition 

Description 

Advance 

ratio 

Load 

factor 

Updated 
(teal gn 
variable 

1 

2 

3 

hover 

forward flight 
maneuver 

0 

0.35 

0.30 

1.0 

1.0 

1.22 


Upper Level Desig n Variables - The design 
model is shown in figure 12. The design variables 


Fig, 11 Lower level of integrated aerody- 
namic-dynamic-structural optimization 
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In equation ( 7 ) , N is the centrifugal force, and EA is 
the cxtensional stiffness. The "coordination" con- 
straints will be defined in detail after the lower level 
optimization is described. 

Lower Level Objective Function - As in- 
dicated previously, a lower level optimization is per- 
formed for each of five cross sections. In each section, 
the objective function is as follows. 


(a) Upper level 



(b) Lower level 


Fig. 12 Design model for integrated aero- 

dynamic/dynamic/structural optimization 

(depicted on the blade schematic) include Elf, El/, 
GJ, EA, and M (at die three spanwise locations 
indicated by the circled numbers), Uirce tuning masses 
and locations, c r , X=cr/ct , yy. and ©tw 

U pper Level Constraints - The upper level 
constraints are the same as diose in the AD procedure 
with the addition of strain constraints and 
"coordination" constraints which assure proper link- 
age and consistency between die upper mid lower lev- 
els. These constraints arc imposed at several span- 
wise locations (currently five) denoted in figure 12(a) 
by the numbered boxes. The locations are die blade 
rooL point of taper initiation, tip, at die two loca- 
tions where the airfoils change (r/R = 0.85 and r/R = 
0 . 95 ). The constraints imposed on die average axial 
strains e x are as follows: 

£ x ! £ a “ 1 - 0 (6) 


2 

F=((m f -Ki;)/Ei;) + 

((EI/ - El* ) / El/ ) + 

((gj-gj*)/gj*)‘ + 

((M-5M')/8M*) 2 (8) 

In equation 8 a starred quantity ( )* denotes a design 
variable from die upper level and 8 represents a 
fraction of the upper level mass M*. The box beam 
designed in the lower level represents structural mass, 
therefore only the structural mass portion of M is 
matched in the lower level. Typically, 8 is chosen to 
be 0.6 which would require the structural mass of the 
box in die lower level to be 60 percent of the total 
upper level mass. The cross sectional properties 1 /, 
I|-, and J arc computed according to engineering theory 

of drill- walled beams 21 from the dimensions of the 
section (b.h.t^.tj, and t 4 ) shown in figure 12(b). 

The overall dimensions b and h are determined from a 
geometrical procedure 22 which fits a rectangle of 
maximum area inside the airfoil section. The airfoil 
section size is determined from the thickness-to-chord 
ratio (t/c) and the local chord (c). The dimension ^ 
represents die local thickness of the airfoil (t^c+t/c). 


l ower Level Desig n Variables - At each 
cross section, die lower level design variables arc die 
four wall thicknesses lj, t2> 13, and t 4 as shown in 
figure 12b for a total of 20 design variables. 

l ower Level Con straints - At a given 
secrion, die lower level constraints are enforced on the 
lower level cxtensional stiffness, cross sectional 
stresses, and wall thicknesses. The extensional 
stiffness constraint requires that the lower level E A 
be greater dian the upper level EA* and is given by 

1 - EA / EA* < 0 (9) 


where £a is the allowable strain and 
e x = N / EA 


( 7 ) 


9 



where EA is the lower level extcnsional stiffness and 
EA is the upper level extensional stiffness at the 
given cross section. It is noted that EA appears in a 
constraint rather than in the objective function 
(equation 8) where the other stiffnesses appear. This 
was done for the following reason. The role of EA in 
the upper level is limited to satisfying the strain 
constraint (equation 7). The lower level is responsible 
only for assuring that the value of EA is at least as 
large as the value needed in the upper level -- close 
matching of EA to EA* is not required. 


,n ^ eam shear str ess analysis. The 
terms 5 and i\ are the principal/centroidal coordinates 
o the appropriate points in the cross section (figure 
13). The moments and the torque,and forces are 
computed m the upper level and passed to the lower 


An additional set of constraints is imposed on the 
lower level design variables to assure that the section 
remains a thin-walled section. These constraints are 


The stress constraints arc evaluated at die comers 
and midsidcs of the cross section shown in figure 13. 



Fig. 13 Locations in cross section for 
stress calculations 

The stress constraints have the following form 

v ( a x.Oxy)/o a -l50 (10) 

where V(ox> a xy) is Von Mises stress measure given 
by 

V(o x .o xy ) = + 3 aly (li) 

where o x , o xy and a a are the axial bending stress, 

shear stress and allowable stress, respectively. The 
stresses a x and c X y are given by 

°x = (M f / If )5 + (M/ / 1/ )q + N / A (12) 


•i / o <().| i = 2and4 (14a) 

tj / h < 0. 1 i = 1 and 3 (14 b) 

where b and h are the width and height of the cross 
section respectively. 




Lyve|s - As mentioned in the section on upper level 
constraints, the coordination between the upper and 
lower levels is implemented by upper level con- 
straints. these constraints are imposed to encourage 
changes in the upper level design variables which 
promote consistency between the stiffnesses from the 
upper and lower levels. Specifically, these 
constraints (one at each of five spanwise stations, 
figure 1 2) have the form 


9 = F-(l + e)F 0 <0 


(15) 


where F is given by equation (8) for the current upper 
level design variables, F 0 is the latest value of the 

lower level objective function, and e is a specified 
tolerance. 

... ik rtvmiv* »f coordination r^ trnim 
1 he derivative of the constraint in equation (15) is re- 
quired as part of the optimization (as are the deriva- 
tives of all the constraints and the objective func- 
tions). The derivative of the coordination constraint 
involves some complexities which are illustrated 
herein. In view of equation (8) a constraint g may 
be written as 6 3 


o w =M T /(2A c l) + 52 


(13) 


9=g{F(X),F 0 (X.Q(X))} 


(16) 


where Mf is the flapwisc moment, M/is die lag 
moment. If is the flapwise bending inertia, 1/ is the 
lag bending inertia, Mj is the torque at the section; 
A c is the area of the cross section enclosed by die 
centerlines of the walls; t is die side wall thickness; 
S is the vertical shear force, Q is the usual first 


where X is the set of upper level design variables and 
K} is die set of behavioral quantities which are 
computed in the upper level and used in the lower 
level ( i.e, Mp M/ , N ). From equation (16) 

ag = _agaF + _dg_a^ a^a Fo aQ 

dx dF ax 3F 0 ax ap 0 ao ax 


l 0 



where = i, if. is obtained by differentiating 

3F 9X 

equation (8>. |t— (H-e) and || >'« 

derivatives of an optimum design with respect to pa- 
rameters 23 and is a behavior sensitivity 

ax 

derivative computed in the upper level. 

Overall Organi z ation of APS ffWttly tf 

The overall system is shown schematically in 
figure 14 and outlined herein. The upper level anal- 



Fie. 14 Overall ADS system implementa- 
tion 

ysis, as shown in figure 10, is executed for the 
current set of design variables Ihis analysis provides 
all of the information needed to calculate the objective 
function and constraints with the exception of die co- 
ordination constraint (equation 15). The lower level 
optimization (figure 11) is then performed for each of 
the five cross sections shown in figure 12 to match 
(as closely as possible) the current upper level 
bending and torsional stiffnesses. The lower level 
optimization also provides the objective function, l"o 
needed for the coordination constraints. These final 
constraints, together with the remaining information 
from the upper level, permit the upper level optimiza- 
tion to take place. This describes one cycle of the 


procedure. Ihc procedure is repeated for additional cy- 
cles until convergence is achieved. 

Concluding Remarks 

This paper has described a joint activity involv- 
ing NASA and Army researchers at the NASA 
Langley Research Center to develop rotorcraft 
optimization procedures which integrate appropriate 
disciplines and account for all of the important 
interactions among the disciplines. The disciplines 
involved include rotor aerodynamics, rotor dynamics, 
rotor structures, airframe dynamics, and acoustics. 
The work is focused on combining the five key 
disciplines listed above in optimization procedures to 
satisfy multidisciplinary design requirements. 
Fundamental to the plan is a three-phased approach. 
In phase I, the disciplines of blade dynamics, blade 
aerodynamics, and blade structure will be closely 
coupled while acoustics and airframe dynamics will be 
decoupled and accounted for as effective constraints on 
llic design for the first three disciplines. In phase 2, 
acoustics is to be integrated with the first three 
disciplines. Finally, in phase 3, airframe dynamics 
will be fully integrated with the other four 
disciplines. The paper described the plan and 
summarized the most recent results and progress in 
the work. 'ITiese include results for optimal placement 
of tuning mass for reduction of blade vibratory shear 
forces, results for integrated aerodynamic/dynamic 
optimization, efforts for validation of the procedures 
and the formulation of a multilevel integrated 
aerodynamic/dynamic/structural optimization pro- 
cedure. The results demonstrate the potential of 
optimization in design of future rotorcraft, both from 
the standpoint of efficiency of the process as well as 
potentially improved products. The results 
demonstrate that there are significant opportunities 
awaiting analytical designers who pursue 
interdisciplinary design approaches. It has been the 
intention of the authors and the Langley team, in 
general, to develop the research with the long-term 
needs of the rotorcraft industry in mind. This paper 
already reflects a significant industry influence as a re- 
sult of modifications following industry visits. It is 
planned to continue this dialogue with the anticipa- 
tion of additional industry interaction during future 
phases of this research. 
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